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On the interpretation of age and chemical composition of composite 
stellar populations determined with line-strength indices. 
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ABSTRACT 

We study the simple-stellar-population-equivalent (SSP-equivalent) age and chemical 
composition measured from the Lick/IDS line-strength indices of composite stellar popu- 
lations (CSP). We build two sets of ~30000 CSP models using stellar populations synthesis 
models, combining an old population and a young population with a range of ages and chemi- 
cal compositions representative of early-type galaxies. We investigate how the SSP-equivalent 
stellar parameters of the CSP's depend on the stellar parameters of the two input populations; 
how they depend on V-band luminosity-weighted stellar parameters; and how SSP-equivalent 
parameters derived from different Balmer-line indices can be used to reveal the presence of a 
young population on top of an old one. We find that the SSP-equivalent age depends primarily 
on the age of the young population and on the mass fraction of the two populations, and that 
the SSP-equivalent chemical composition depends mainly on the chemical composition of 
the old population. Furthermore, while the SSP-equivalent chemical composition tracks quite 
closely the V-band luminosity weighted one, the SSP-equivalent age does not and is strongly 
biased towards the age of the young population. In this bias the age of the young population 
and the mass fraction between old and young population are degenerate. Finally, assuming 
typical error bars, we find that a discrepancy between the SSP-equivalent parameters deter- 
mined with different Balmer-line indices can reveal the presence of a young stellar population 
on top of an old one as long as the age of the young population is less than ~2.5 Gyr and the 
mass fraction of young to old population is between 1% and 10%. Such disrepancy is larger 
at supersolar metallicities. 
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1 INTRODUCTION 

The knowledge of age and chemical composition of stars in early- 
type galaxies is a fundamental piece in the puzzle of galaxy forma- 
tion and evolution. For a long time, optical-wavelength studies in 
this direction have been hampered by the age-metallicity degener- 
acy: an age variation of a factor of ~3 mimics a metallicity varia- 
tion of a factor of ~2 in the spectra of old stellar populations (e.g., 
Faber 1973; O'Connell 1986; Worthey 1994). The effort of vari- 
ous authors during the past two decades culminated in the work of 
Worthey (1994), who showed that age and metallicity can be disen- 
tangled by the joint use of pairs of line-strength indices, one metal- 
line and one Balmer-line index, measured from the optical spec- 
tra of galaxies. A system of line-strength indices was defined (the 
Lick/IDS system, see Burstein et al. 1984; Worthey et al. 1994 and 
references therein; Worthey & Ottaviani 1997) and is now widely 
used in order to determine the age t, metallicity [Z/H] and abun- 
dance ratio [E/Fe] of stars in galaxies ([E/Fe] is defined in Trager 
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et al. 2000a as a way of parameterising deviations from the solar 
abundance pattern). 

In practise, one compares the indices measured from the opti- 
cal spectrum of a galaxy to their values predicted by stellar popula- 
tions models (provided for example by Worthey 1994; Vazdekis 
1999; Bruzual & Chariot 2003; Thomas, Maraston & Bender 
2003). The stellar t, [Z/H] and [E/Fe] of the galaxy are the ones 
of the model whose indices best agree with the measured ones. 
Because each model consists of a single-burst stellar population 
(SSP) whose stars, unlike in real galaxies, all have the same t, 
[Z/H] and [E/Fe], the derived stellar parameters are labelled as SSP- 
equivalent. We will refer to them as r SSP , [Z/H] SSP and [E/Fe] SSP . 

Early-type galaxies are the most massive stellar systems for 
which the SSP approximation seems to hold. Therefore, many au- 
thors have measured their line-strength indices in order to deter- 
mine their stellar content (see for example Trager et al. 2000b; 
Caldwell, Rose & Concannon 2003; Denicolo et al. 2005; Thomas 
et al. 2005; Clemens et al. 2006). However, many results suggest 
that recent star formation occurred in these galaxies (Trager et al. 
2000b; Yi et al. 2005), so that a small fraction of their current stellar 
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Table 1. Input parameters for the composite stellar populations 



parameter values 

t\ (Gyr) 10.0,13.0 

[Z/Hh -0.15,0.0,0.15,0.3 

[E/Fe]i 0.0,0.15,0.3,0.45 

f 2 (Gyr) 1,1.4,1.8,2.5,3.4,4.7 

[Z/Hh 0.0,0.15,0.3,0.45 

[E/Feh -0.15,0,0.15,0.3 

M =M 2 /Mi 0.0,0.001 ,0.005,0.01 ,0.025,0.05,0. 1 ,0.2,0.35,0.5 



Stellar parameters are labelled as "1" and "2" corresponding to the old pop- 
ulation (SSPi ) and the young population (SSP2) respectively. The last row 
lists the values adopted for the mass fraction between the two populations. 

mass formed few Gyr ago. It is natural to wonder about the meaning 
of the line-strength indices analysis under such circumstances (i.e., 
in the presence of more than one SSP in the same galaxy). How do 
the SSP-equivalent parameters relate to the average properties of a 
galaxy and to the ones of the many SSP's that it hosts? And what 
do we actually learn from SSP-equivalent parameters? 

Previous authors have already looked into this problem. Using 
a limited number of models of composite stellar populations (CSP), 
Trager et al. (2000b) found that the SSP-equivalent age is heavily 
biased towards the age of the young stars present in a galaxy. In 
this paper we address the same questions in a more systematic and 
extensive way from the point of view of the stellar population mod- 
els. 

We build two sets of CSP's by using the models of Bruzual 
& Chariot (2003, hereafter BC03) and Worthey (1994, hereafter 
W94). Each dataset contains ~3xl0 4 CSP models composed of 
one old and one young SSP. Different CSP's correspond to dif- 
ferent stellar parameters of the parent SSP's. The old SSP (SSPi) 
is always chosen to be more massive than the young one (SSP 2 ), 
as inferred in many early-type galaxies (e.g., Trager et al. 2000b; 
Leonardi & Worthey 2000; Jeong et al. 2006). We analyse the line- 
strength indices of the CSP's and derive the corresponding f SSP , 
[Z/H]ssp and [E/Fe]ssp as is usually done for observed galaxies. 
We then analyse how the result depends on the input parameters 
(ri, [Z/H],, [E/Fe],, t 2 , [Z/H] 2 ,[E/Fe] 2 and pL=M 2 /M x where M is 
the stellar mass) and on the luminosity-weighted properties of the 
CSP's. We restrict our study to systems composed of only two 
SSP's because this case is still relatively easy to treat and might 
be a reasonable first-order approximation for early-type galaxies. 
In Sect|5|we explain in some details the construction of the two 
datasets, we present the results in Sect[3] and finally draw some 
conclusions. 



2 MODELS OF COMPOSITE STELLAR POPULATIONS 

We define a set of old and young SSP's in TableQ All the possi- 
ble combinations of one old and one young SSP for each of the 
values of the mass fraction /j give the 30720 CSP models that 
form each of the datasets. A dataset consists of a 7-dimensional 
grid in the parameters (t u [Z/H],, [E/Feh, h, [Z/H] 2 , [E/Fe] 2 , 
/i). At each position in the grid we store a vector that con- 
tains the line-strength indices of the two input SSP's and of the 
resulting CSP, the luminosity-weighted stellar parameters (?lw. 
[Z/H] LW , [E/Fe] LW ) and the SSP-equivalent stellar parameters (t SSP , 
[Z/H] SS p,[E/Fe]ssp). 




[MgFe] 




[MgFe] csp 

Figure 1. Top panel: BC03 parent SSP's of solar [E/Fe] plotted on top of 
BC03 model grid; filled triangles represent the young SSP's (SSP 2 ); empty 
circles represent the old SSP's (SSPi) which dominate the mass of the 
CSP's (see Table [TJ- Bottom panel: distribution of the BC03 CSP's built 
from the solar-[E/Fe] SSP's. 



The luminosity-weighted stellar parameters are defined in 
terms of the V-luminosity of the parent populations at the respec- 
tive ages (i.e., taking into account the change in mass-to-light ratio 
as a population ages). 

The SSP-equivalent parameters are determined by comparing 
the Lick/IDS line-strength indices Mg/?, Fe5270 and Fe5335 and 
a Balmer-line index of the CSP to their values according to mod- 
els (BC03 or W94 depending on the dataset). As Balmer-line in- 
dex we use alternatively HyS, Hy A , Hy F , K6 A or K6 F , obtaining 
SSP-equivalent parameters for each of them separately. These will 
then be labelled according to the Balmer-line index used in deriv- 
ing them (e.g., t H p is the SSP-equivalent age derived using Mgfo, 
Fe5270 and Fe5335 and H/3). The use of different Balmer lines 
is important because they respond differently to the presence of 
a young stellar component on top of the old one (Schiavon et al. 
2004). 

Although W94 and BC03 models are available only with 
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Table 2. Covariance coefficients between SSP-equivalent parameters and input parameters 



logf! logf 2 [Z/H], [Z/H] 2 [E/Fe]i [E/Fe] 2 /j 



logr S sp 0.05,0.10 0.51,0.56 -0.14,-0.07 0.09,0.12 -0.01,0.01 -0.01,0.00 -0.57,-0.(55 
[Z/H]ssp 0.02 , 0.06 -0.19,-0.32 0.78 . 0.85 0.09 , 0.16 -0.01,0.01 0.00 , 0.01 0.25 
[E/Fe] SS p 0.00,0.01 -0.11,-0.01 0.02.0.07 -0.06,-0.04 0.85,0.87 0.23 -0.19,-0.17 



Each entry is the range within which the covariance coefficient varies when changing Balmer-line index (or the value of the coefficient if this does not vary). 

The largest coefficients for each SSP-equivalent parameter are given in italics. 



[E/Fe]=0, TableQcontains also parent SSP's of non-solar [E/Fe]. In 
these cases we correct the line-strength indices given by the models 
according to the [E/Fe] variations. We then use the corrected val- 
ues when both building the CSP models and measuring their SSP- 
equivalent parameters. The correction scheme is the one described 
in Trager et al. (2000a) but improved by the use of new response 
functions computed and kindly provided by G. Worthey. For details 
see Trager, Faber & Dressier (2006). 

FigQ shows the distribution of the solar-[E/Fe] parent SSP's 
drawn from the BC03 models and of the resulting CSP's on the 
plane [[MgFe],H/J], where age and metallicity are efficiently de- 
coupled. The points are plotted on top of the BC03 solar-[E/Fe] 
model grid. The CSP models cover most of the area where early- 
type galaxies have been observed to lie (e.g., Trager et al. 2000a; 
Denicolo et al. 2005; Thomas et al. 2005). 



3 RESULTS 

We calculate the covariance coefficients between the SSP- 
equivalent stellar parameters and the input parameters in order to 
understand which input parameters are mostly driving the varia- 
tions in fssp, K/H]ssp an d [E/Fe]ssp- Table|2|shows the result of this 
calculation performed on the BC03 dataset. The covariance coeffi- 
cients do not change much when changing Balmer-line index. In the 
table we give the range within which each coefficient varies when 
changing Balmer-line index. Furthermore, we have verified that the 
W94 dataset gives the same result. The following comments apply 
therefore to all Balmer-line indices and to both datasets. 

• The variations in fssp are mostly driven by variations in f 2 , the 
age of the young population, and fj, the mass fraction, while other 
parameters play a secondary role. The sign and absolute value of 
these two covariance coefficients clearly show the strong degener- 
acy between f 2 and ji: the same /ssp can result from a small mass of 
young stars or a sufficiently large mass of older stars. 

• The variations in [Z/H]ssp are by far dominated by variations 
in [Z/H]i, the metallicity of the old population. The mass fraction 
/i and the age of the young population f 2 also play a relevant role 
with the usual degeneracy. The latter correlations must be (at least 
partially) due to the fact that in the dataset [Z/H] 2 is on average 
larger than [Z/H] i . However, we have verified that the covariance 
coefficients between [Z/H]ssp and \i and between [Z/H]ssp and f 2 
remain significantly larger than zero when considering a subset of 
models with [Z/H][ and [Z/H] 2 sampled in identical ways (in par- 
ticular, the covariance coefficients drop by a factor of ~2 and ~1.3 
respectively). 

• [E/Fe] ssp seems to be varying mostly because of variations in 
the abundance ratios of the two parent populations, with the old, 
massive population being dominant. As for [Z/H]ssp, the correla- 



tion with yti and f 2 is only in part due to the different sampling of 
[E/Fe]! and [E/Fe] 2 . Using a subset of models with identical sam- 
pling of [E/Fe]i and [E/Fe] 2 reduces the covariance coefficient with 
H and increases the one with t% by a factor of ~3. 

Covariance coefficients highlight which of the input parame- 
ters play the dominant role in determining the variation of the SSP- 
equivalent ones. It is also interesting to see how the latter relate to 
the average properties of the model CSP's. Fig|2| shows the com- 
parison between the H/3-based SSP-equivalent parameters and the 
V-band luminosity-weighted ones. The behaviour is substantially 
the same when using different Balmer-line indices. 

[Z/H]ssp and [E/Fe]ssp seem to track quite closely [Z/H]lw 
and [E/Fe] LW respectively; strong deviations are observed only for 
the youngest f 2 . On the other hand fssp is always much smaller than 
?lw and lies somewhere between the latter and f 2 . This effect was 
already known (Trager et al. 2000b). Its explanation is that the de- 
termination of /ssp relies primarily on Balmer-line indices (see the 
model grid in Fig0. These are dominated by young stars and there- 
fore f S sp is strongly biased towards the age of the young stellar 
component. As Fig[2|illustrates, the younger SSP 2 the stronger this 
bias is. 

Fig|5| demonstrates that it is not correct to use f S sp as an esti- 
mate of when a galaxy formed its stars (yet, this is often done; see 
for example Clemens et al. 2006). A fair statement would be that 
fssp is a Balmer-line-weighted age and it should always be kept 
in mind that such age is strongly biased towards the age of young 
stellar components. Furthermore, as highlighted by the covariance 
coefficient and confirmed by Fig|5| the effect of f 2 and /j on f S sp 
is degenerate. An increasingly older SSP 2 can produce the same 
fssp as long as fj is properly increased (in Fig|2|p increases towards 
decreasing f S sp and /lw)- 

As mentioned, SSP-equivalent parameters derived from differ- 
ent Balmer-line indices behave substantially in the same way (i.e., 
FiglU looks roughly the same for all of them). However, different 
Balmer-line indices are sensitive to the presence of young stars at 
different levels (Schiavon et al. 2004). Because of this fssp, [Z/H]ssp 
and [E/Fe]ssp (of a CSP) computed with different Balmer-line in- 
dices will not be in agreement. Fig[3] illustrates this concept for a 
subset of the CSP models chosen to have solar chemical composi- 
tion and fi = 13 Gyr. It can also be seen that the difference between 
H/3- and Hy^-based SSP-equivalent parameters goes to zero for fi 
approaching and 1 and peaks between these two extremes at a 
position dependent on f 2 . Furthermore, Fig|5|shows once more the 
degeneracy between f 2 and /i. The same difference between, for ex- 
ample, t H p and t HyA can be caused by increasingly older SSP 2 's as 
long as the mass fraction yi is sufficiently increased. 

For clarity Fig[3] shows only a subset of the CSP models. A 
similar trend is anyway observed in the whole sample (and in W94 
dataset), showing that it is possible to detect the presence of a 
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Figure 2. H/3-based fssp, [Z/H]ssp and [E/Fe]ssp plotted versus the respec- 
tive V-band luminosity-weighted quantities using BC03 dataset. The dashed 
line of each plot is the identity line. The colour codes ti . 



young stellar component on the basis of the disagreement between 
SSP-equivalent parameters obtained with different Balmer-line in- 
dices. We actually expect that more dramatic disagreements in, for 
example, rssp are accompanied by larger differences in [Z/H]ssp and 
[E/Fe] SSP . This is indeed observed and showed in Fig[3] where the 
difference between H/3- and Hy A -based [Z/H] S sp and [E/Fe] SSP is 
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Figure 3. Difference between H/?- and Hy^ -based SSP-equivalent parame- 
ters as a function of //. The plots are obtained using the BC03 dataset and 
solar [Z/H] and [E/Fe] for both SSPi and SSP 2 . The age of SSP! is fixed at 
13 Gyr. 



plotted versus the difference in f SSP . In particular the age-metallicity 
plot shows a very tight relation. This could be used in order to test 
the correctness of one's results. 

Different B aimer lines can be used as a diagnostic for the 
presence of a young stellar component only as long as the differ- 
ences in the SSP-equivalent parameters are larger than the obser- 
vational errors. These are typically of 0.1 on the logarithm of t SSP 
and on [Z/H] SSP and of 0.05 on [E/Fe] SSP (see for example Trager 
et al. 2000a; Thomas et al. 2005). Fig|4]shows that with these er- 
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Figure 4. Difference between H/?- and Hy^ -based [Z/H]ssp (top) and 
[E/Fe]ssp (bottom) plotted versus the difference in fssp- Each point cor- 
responds to a BC03 CSP model. The presence of a young stellar population 
on top of an old one causes SSP-equivalent parameters based on different 
Balmer-line indices to disagree. This effect must be and indeed is observ- 
able simultaneously in fssp. [Z/H]ssp and [E/Fe]ssp- In particular, points 
seem to be distributed along a very tight relation in the age-metallicity 
plane. 



rors fssp measurements are the most efficient in revealing a young 
component, allowing the detection of SSPa's younger than ~2.5 
Gyr. As suggested by Fig|3| this is however possible only within 
a certain range of yi. The actual range depends on tj but we find 
it to be roughly between 1% and 10%. [Z/H] SSP and [E/Fe] SSP are 
only sensitive to SSP2's younger than ~1.5 Gyr with /i between 2% 
and 10%. It is important to stress that for fi >10% there is no de- 
tectable difference between the SSP-equivalent parameters derived 
from different Balmer-line indices. At these values of fi the Balmer- 
line indices are so heavily dominated by the younger populations 
that they all "see" the same age, which is very close to the age of 
the young population. 

Fig[3]and Fig|3]show another interesting feature: the largest 
difference between fH/3 and tn 7A in Fig|5|is small compared to the 
one in Fig[3] Recall that Figgis relative to CSP's where both SSPi 
and SSP2 have solar chemical composition, while Fig^Jrepresents 
the whole BC03 sample, with [Z/H] growing significantly above 
solar. The two figures suggest that [Z/H] plays an important role 
with respect to the difference t- l y } -t }lyA . Fig|5| shows that indeed 
high (and therefore more easily detectable) differences in fssp oc- 
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curr only at high metallicities. Similar plots hold for the difference 
in [Z/H] and [E/Fe]. We therefore find that the use of more than one 
Balmer-line index can reveal the presence of a young stellar popu- 
lations, but that this is possible only for a small range of fi and \i 
and depends also on the metallicity of the populations. 

We would like to remind the reader that a disagreement be- 
tween Balmer-line-based SSP-equivalent parameters, in principle 
revealing the presence of a young stellar component, could also re- 
sult from the approach used when analysing the data. In particular, 
it is important to remember that different Balmer-line indices re- 
spond differently to variations in [E/Fe]. Thomas et al. (2004) and 
Thomas & Davies (2006) pointed out that this causes a discrep- 
ancy between SSP-equivalent parameters determined from differ- 
ent Balmer-line indices when using as a comparison models with 
solar [E/Fe] only. This effect could mimic the presence of a young 
stellar component. However, no such problem should occur when 
using models that account properly for [E/Fe] variations, as was 
done here. 

Another delicate point when using Balmer-line indices is their 
increase caused by hot star populations like blue horizontal branch 
stars or blue stragglers (Maraston & Thomas 2000; Trager et al. 
2005). Although this is not an issue for the present study, where 
we do not explore the metal-poor regime at which these stars are 
expected to be found, this problem should always be kept into con- 
sideration when dealing with real data. 



4 CONCLUSIONS 

We have built two datasets of composite stellar populations (CSP) 
using Bruzual & Chariot (2003) and Worthey (1994) models. Each 
CSP model in the datasets consists of an old single-burst stellar 
population (SSPi) and a younger, less massive one (SSP2). We 
have investigated how the SSP-equivalent parameters determined 
by measuring the Lick/IDS line-strength indices of the CSP's de- 
pend on the stellar parameters of SSPi an d SSP2. By means of co- 
variance coefficients we have found that, regardless of the particu- 
lar stellar populations models used in building the CSP's and of the 
Balmer-line index used for the analysis: f SSP , the SSP-equivalent 
age, depends primarily on t%, the age of the young population, and 
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p., the mass fraction between the two populations; variations in 
[Z/H] S sp, the SSP-equivalent metallicity, are mostly driven by vari- 
ations in [Z/H]i, the metallicity of the old population; and [E/Fe]ssp 
the SSP-equivalent abundance ratio, depends mainly on [E/Fe] [ , the 
abundance ratio of the old population. 

Furthermore, we have found that [Z/H] SSP and [E/Fe] SSP track 
quite closely the V-band luminosity-weighted metallicity and abun- 
dance ratio ([Z/H]lw and [E/Fe]L\v) except in case of very young 
(and significantly massive) SSP2. On the other hand, / SSP does not 
follow /lw. being strongly biased towards the the age ti of the 
young population. The SSP-equivalent age /ssp is simply a Balmer- 
line-weighted age and should not be interpreted as the time passed 
since the formation of most of the stars in a galaxy. 

Finally, as found by Schiavon et al. (2004), using more than 
one Balmer-line index can reveal the presence of a young stellar 
component on top of an old one. In this case, SSP-equivalent pa- 
rameters derived from different Balmer-line index give discrepant 
results. This is true however only for values of p between 1% and 
10% and for t 2 <2.5 Gyr assuming typical errors on f S sp, [Z/H] SSP 
and [E/Fe]ssp- Furthermore, these discrepancies are higher at su- 
persolar [Z/H] SSP . Finally, this method does not appear to break the 
degeneracy between the age and the mass fraction of the young 
population, especially when considering the size of the typical er- 
ror bars. In this respect, what is really needed is an age-sensitive 
index dependent on the age of the old stellar population (i.e., RGB 
stars), to be used in combination with Balmer-line indices. 
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